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ABSTRACT 

A time series of pressure variations was obtained by digitizing 
analogue records of bottom pressure, recorded at Pt. Sur, California, 
during the periods November, December, and January, 1947-1948, and 
June, July, and August, 1949. The power spectrum of the pressure vari- 
oe was found via the amplitude spectra obtained by means of the 
fast Fourier transform. This spectrum was then smoothed by averaging 
over several frequency bands and confidence limits were assigned. 

Surface wave power spectra were inferred from the pressure power 
spectra and significant wave heights were estimated from these, after 
applying a procedure designed to remove noise. 

Values of significant height on the bottom and at the surface 
are displayed in tabular form and plotted versus time for the selected 
periods January, 1948, and August, 1949. The spectral patterns of wave 


activity show the transient nature of the energy due to storm swell 


and local winds. 
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TABLE OF SYMBOLS AND ABBREVIATIONS 


Definition 


the maximum frequency about which the 


analysis can gain information 
time 


sampling interval 


L/ Be. elementary spectral band width 


energy 
pressure amplitude 


number of observations 


IBM computer subprogram, a fast Fourier 


transform 

complex spectral amplitude 
frequency of Fourier amplitude 
time of the pressure data point 
Energy in band width Aft 
arbitrary number of band widths 


average spectral amplitude in band 
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Symbol 


RO 


H3B 


H3S 


PCT 


var 


Definition 


summation of the fA y for the ar 


position in time 


radius of the ordinate on the wave 
record 


measure of arc length 
significant height of pressure amplitudes 
measured at the recorder (in this case on 
the bottom) | 


significant height of the surface waves 
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variance 


degrees of freedom 
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I. Objectives 


The randomness of the sea has puzzled investigators for centuries. 
Before forecasts of the behavior of the sea can be made, there must 
exist a more complete knowledge of the many-faceted energy relationships 
between the air and the sea. Most modern authorities feel that the 
true behavior of the sea can be best approximated by considering the 
behavior of each of its spectral components. , 

In addition, with the increasing use of major structures situated 
on the sea floor, it becomes more important to know about the possibil- 
ity of pulsating water motions, with particular frequencies, and changes 
in the magnitude of the pulsations. The knowledge of bottom pressure 
activity is of extreme importance in mine warfare. In view of these 
and other needs, this paper has the following objectives: 

A. To display the spectral patterns of energy density of waves 

and associated bottom pressures versus time for the analogue 
wave records recorded at Pt. Sur, California, during the 
periods: 

1. November, December, and January, 1947-48, (winter), and, 

2. June, July, and August, 1949, (summer). 

The spectra are to be obtained by digitizing the analogue data 
and performing a spectrum analysis on the resulting time series. 

B. To compare the results from the summer and winter periods; 

C. To provide background data for other studies, such as; 

1. Consideration of effects of shoaling and refraction in 
obtaining a deep water spectrum; 

2. Wave energy attenuation studies; 

3. Locating of sources of wave activity reaching the Pt. Sur 


area; and, 


ll 


4. Examining the analysis for peculiarities such as favored 
frequencies that might indicate the influence of local 


topography. 
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II. Related Work 


Munk, Snodgrass, and Miller (1) displayed contours of equal power 
density, E(f,t) on a frequency time plot and used the spectral peaks' 
arrival times to estimate distances to the areas of wave generation. 
They confirmed the relationship by determining direction and distance 
in a separate analysis. 

The use of a digitizer to convert the analogue data of the wave 
records is discussed by Moskowitz, Pierson, and Mehr (2). Their 
records were part of the data used in "the problem of developing 
numerical wave forecasting procedures for the North Atlantic Ocean", 


now functioning as a part of ASWEPS. 
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III. Data Used in This Study 


The wave pressure records used for this study were recorded by a 
University of California Mark III bottom pressure recorder, anchored 
on the bottom in sixty-eight feet of water, measured from MLLW. A 
comprehensive discussion of bottom pressure recorders and the manner 
in which they act as a low pass filter was presented by Munk (3). 

The major advantage of choosing these somewhat dated records is that 
they provide a long and nearly continuous history from which patterns 
of wave activity may be shown. Johnson (4) has also published the re- 


sults of another kind of analysis of these same records. 
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IV. Spectral Analysis 


A time series of pressure variations is analyzed using the fast 
Fourier transform, a technique developed by Cooley and Tukey (5), in- 
corporated in the IBM computer subroutine HARM to obtain the amplitude 
spectrum. The Fourier amplitudes are then used to find the power spectrum 
which is smoothed by averaging over several frequency bands. The entire 
procedure and several ee processes have been incorporated into a 
computer program written in FORTRAN IV. 

A basic step in the analysis is selection of a suitable sample inter- 
val. For this study the primary interest was in wind waves; and, both 
theory and observation indicate that, due to filtering, very little 
activity at depths of sixty-eight feet is found in waves of frequencies 
greater than 0.25 a The appropriate sampling interval for this 


frequency is A t = 2 seconds, given by 


ec 
‘max 2DAt (1) 


The fast Fourier transform requires 2” (n a positive integer) 
discrete data points for each analysis. With a 2-second sampling inter- 
val, the largest n possible was 9 (512 data points) equivalent to 17 
minutes and 4 seconds. 

As a check on the spectral analysis, an independent calculation 
of the total energy was made directly from the pressure amplitudes, by 


N 2 
E = 2 var(A,) = 2/N=, AL (2) 


where, A. = pressure amplitude, and, 


M = 512 is the number of observations of A. 
The Fourier series analysis, HARM, produces a complex amplitude 


spectrum which represents the points of the time series by the sum of 


the same number of discrete frequencies. The complex amplitudes are 


given by 


51 
eee om, (3) 
k= 


IN e @ 
where, denotes Fourier transformation, 


rh 
il 


j Af is frequency, and, 


ct 
il 


- k At is the time of the pressure data points. 
The energy in each band of width Af is proportional to 


Bee Eee A. 7 (4) 
J J J 








A 
where the factor of two enters because of the symmetry of oF about 


j = 256, and the total energy is proportional to 


“ 
JI 


6 
a,| 7 2 me (5) 


S12 | 
1 2 


Cute 
il 


Because the sea surface disturbance is a random variable, the 
spectral values ACE.) are also random variables. Thus equation (3) 
holds only for averaging over many realizations. The energy is also 
a random variable and equations (4) and (5) are only estimates of the 
respective energies. The confidence limits for these estimates can 
be determined in the manner presented by Pierson and Marks (6) and 
discussed here below. 

The energy in each elementary band width has been determined from 
only two independent bits of information, and the estimates are distri- 
buted according to the Nos distribution with two degrees of freedom. 
By averaging a number, say m, of elementary bands of width A f (see 
figure 1), one obtains an estimate of the energy which is more stable 


than the estimates of the contributing narrow bands. 
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Figure 1. Bi: representing ME in band width m Af 


These estimates have a a distribution with 2m degrees of freedom. 


Confidence limits may be assigned utilizing this distribution of B, 


Band averaging was done over 10 and 16 By The total energy per band 
of width m A f was computed by 
=l+m/2 
AE(£,) =  |1/m F m= 25, (6) 
j=l-m/2+1 


where, fy is considered the center frequency of the band, and, 


m * number of bands. 


The cumulative energy was obtained by summing over m bands, 


M ~ 7/2 M=- "Yo, 
E = DARE) = DB, . (7) 
2 VG 


The average energy density in the pressure variation in the band 


width mA £ was defined as 
S(£,) = AE(£,)/m Af. (8) 
The average energy density at the surface was then determined by 
S(£,) = 8(£,)/K° (9) 
s* 1 1 s 


where K is the well known pressure response factor, 


1 
m cosh 2 77 d/L ’ ey 


Li 


where, d = depth of water (meter assumed on sea floor), 
and, 


L = wave length at the depth d. 
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V. Performing the Analysis 


The wave records used in this project were analogue records of 
pressure variations recorded on Esterline Angus Co., Inc. chart 4313, 
shown in figure 2, which has curvilinear co-ordinates (t,RQ@). The 
records were digitized using the Calma Co. model 480 digitizer, which 
recorded incremental pairs, in the order, Ax and Ay (of the cartesian 
co-ordinates of the pressure curves), every .01 inch of stylus travel. 
This resulted in the wave form being approximated by a step - function, 
as shown in figure 3, with a secondary wave form including very high 
frequency energy being superimposed on the iecorded wave form. This 


is a source of error and will be discussed further in section VI. 
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Figure 3. Steplike approximation to the recorded trace. 

The digitizer recorded incremental values and it was therefore 
necessary to form the cartesian co-ordinates by summation within the 
computer. These cartesian co-ordinates were then converted to the ap- 
proximately correct t,R@ co-ordinates by employing geometric relation- 
ships. 

Computation of the correct value of t. (the time axis) was critical 
to the conversion. On the time scale, 0.1 inch of record was equivalent 
to the sample interval of 2-seconds. The summation of Ax, x.» repre- 


senting the horizontal travel, was converted to the t,R@ co-ordinates 
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by, 


2 2 
= 4 - = 
t. x, R R ve (11) 


where, x, is the time scale in equally spaced increments of 
At on the cartesian scale, 
y. is the ordinate on the cartesian scale, 
t. is the time in t,RO@ co-ordinates, and, 
R is the constant radius for arcs of constant t,. 


i 


These are shown in figure 4. 


Yao 


Y, 
5 
: “ 
in : 1 


Xi : . 
| 
a R* -Yi 
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Q. = tan y,/ Rice: (12) 


Figure 4. Geometric relations for conversion to t,R®@ 
co-ordinates 


Because 2-second intervals read by the digitizer on the rectangular 
co-ordinates do not in general correspond to the same interval in curvi- 
linear co-ordinates, the following scheme was used to approximate even 
time steps. The cartesian amplitudes used to form the time series were 
extracted from the record when the value of t. exceeded or closely ap- 
proximated a separately incremented test value which represented the 


true position in time; i.e., 
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tert = test or if test -t .001, then 


Yi YG 
where, i is the position in the record, and, 
j is the position in the time series. 
The cartesian amplitudes were then converted to arc length measured 


in scale units; by 


y = RO (13) 


The value of R used in these calculations was determined by graphi- 
cal construction. By repeated measurement, with the digitizer, the 
radius was determined to be 4.36 + .02 inches or approximately 490 scale 


units. 
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VL. Discussion of Results 


One feature of the surface spectra was the characteristic "rapid 
increase in the high frequencies due to the presence of noise .. ." 
as described by Moskowitz, et. al. (2). As indicated in section V, 
the digitizing procedure superimposed a very high frequency waveform 
on the recorded trace, the energy being aliased eventually into the 
frequencies used in the analysis. 

An effort was made to estimate the contribution of this to the 
total energy. The maximum reading error of the digitizer was .01l 
inch on the recording or .0167 feet in the pressure amplitude. As- 
suming that the errors have a uniform distribution (where each event 
is equally likely) between zero and 0.0167 and that the errors in ampli- 
tude are not correlated with the amplitudes themselves, the contribution 
of noise due to digitizing appears to be less than one per cent, usually 
much less. However, even a small amount of noise is critical at the 
higher frequencies, since the energy there is also very small. Moreover, 
there was considerable additional noise not associated with digitizing 
as spectral results of non-digitized data also displayed the rise in 
high frequencies. 

In the absence of more precise information about the total noise 
contribution, it was assumed that most of the energy appearing in the 
spectra at frequencies greater than 0.18 sc. was noise. Much of 
this was then arbitrarily removed in a manner analogous to that employed 
by Moskowitz et. al. (2), wherein it was assumed that after aliasing 
and smoothing over a number of spectral bands, the noise was white. The 
resultant determination of noise was subtracted from all spectral bands 


thus giving spectra with reasonable energy in the higher frequencies. 


23 


All values above 0.25 aneett indicated on computer printouts are beyond 
the range of the sampling interval and are spurious. 

The total energy computed by the sum of the Fourier amplitudes 
agreed well with that found from the variance of the pressure amplitudes. 
The results of band averaging over 16 and 10 elementary band widths were 
essentially the same (in part due to the smoothing to eliminate noise); 
however, the greater detail expected from averaging over fewer bands was 
not entirely lost as can be seen from the selected samples in Appendices 
I and II. 

The corrected values of bottom pressure spectra and inferred surface 
spectra were used to calculate the respective significant heights. The 
significant heights at the surface, for the month of January, 1948, are 
compared with values obtained by Johnson (4) in figure 5. Figures 6 
and / show values of significant height on the bottom and at the surface 
plotted versus time for the selected periods January, 1948, and August, 
1949. A tabular display of significant heights for all periods of this 
investigation are presented in Appendix III. 

An attempt was made to determine distance to the generating area of 
the swell by observing the time lags between spectral peaks of increas- 
ing frequency, as done by Munk et. al. (1). Figures 8 and 9 are contours 
of equal power density, drawn for the characteristic of the logarithm 
of the spectral component of energy density (N = 32) at the surface. 

The very steep slopes are indicative of a relatively close storm area. 
This would be expected for the winter months. The summer months show 
no significant patterns either, perhaps due to the smoothing. In Munk's 
investigation a much more sensitive pressure recorder was used (vibra- 


tron developed by F. E. Snodgrass). The great amount of correlated 
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detail upon which Munk's hypothesis was based occurred between 0.40 
and 0.80 gaara? and therefore is beyond the range of this analysis. 
The ambiguity might be cleared up by more frequent analyses, say 12 


or 6 hours apart vice 24. 
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VII. Conclusions and Recommendations 

Intensity peaks in the summer and the winter occur about a week 
apart. The surface values of significant height in winter are about 
twice those of the summer period. The bottom ratios are much smaller. 
This difference is due to the different distribution of the energy with 
respect to frequency between winter and summer. There seems to be a 
preferred band ae about 0.08 ae for the winter. No 
Similar characteristic appears for the summer months. 

Further investigations to determine distance to the generating 
area should be conducted with a more sensitive digitizing instrument 
located at a shallower depth with analyses more often than daily. These 


procedures would tend to eliminate the high frequency noise. 
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APPENDIX [I 


TABULAR DISPLAY OF SPECTRAL COMPONENTS FOR SELECTED INTERVALS 
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APPENDIX II 


GRAPHIC DISPLAY OF SPECTRAL COMPONENTS FOR SELECTED INTERVALS 
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SUMMARY OF H3B AND H3S FOR NOVEMBER, 1947 


Date Time (PCT) H3B (ft.) H3S (ft.) 
2 1138 - 1156 1.299 42.528 
5 Vie, - 1142 0.436 9.969 
6 Wey - 1205 0.621 7.408 
7 1143 - 1201 0.413 5.008 
8 P23 - 1200 417 2.766 
9 1142 - 1206 1.075 6.531 
10 1150 - 1208 0.995 6.224 
11 1145 - 1203 0.575 2.687 
Ww Wig = 1200 0.754 8.125 

13 1202. -1210 0.597 4.922 
14 1800 - 1818 0.853 7.186 
15 1202 291220 1.022 5.042 
16 1153 - 1211 0.971 8.036 
17 1805 - 1823 0.820 5.255 
18 0854 - 0912 He ileal 10.150 
19 2103-2 0.768 13.533 

20 1459 - 1517 1.148 6.693 

21 1456 - 1514 1.382 5.597 

22 1519 - 1537 1.036 4.355 

23 1617 - 1635 027.27 1.738 

24 1605 - 1623 0.588 2.917 
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Time 


1323 
1340 
1305 
1820 
1804 


1102 
1107 
2009 
0804 
1421 
1418 
1421 
2022 
1420 
1410 


(PCT) 


- 1341 
- 1349* 
- 1323 
- 1838 
- 1822 


- ere 
- 1116* 
- 2027 
- 0822 
- 1439 
- 1436 
- 1439 
- 2040 
- 1438 
- 1428 


* - Record doubled 
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H3B (ft.) 


rp 


Orr Or © kK DN ND NO 


SUMMARY OF H3B AND H3S FOR DECEMBER, 1947 


.718 
-105 
.326 


827 


. ae 


.306 
. Oza 
491 


lw 


194 
. 100 
875 
- 920 
. ie 
-473 


H3S..@re™ ) 


Oi 
.210 
ae 
.085 
Wo 


UW DOW kN 


.341 
2 
.419 
BLS, 
. 902 
226 
~ 482 
10.276 
3.2035 
1.679 


— — 
We Pn Oo fw 


SUMMARY OF H3B AND H3S FOR JANUARY, 1948 


Date Time (PCT) H3B (ft.) H3S (Gite) 
i 0818 - 0836 0.278 1.206 
2 0814 - 0832 1.159 6.309 
3 1428 - 1446 3.28 11.5973 
4, 1422 - 1440 2.795 13 484 
5 0820 - 0838 1.840 5.140 
6 1432 - 1450 0.951 4.811 
7 1610 - 1628 1. 572 7.384 
8 1005 - 1023 1.804 6.387 
9 0943 - 1001 2.458 9.805 

10 1634 - 1652 1.492 5.567 

11 1610 - 1628 1.274 5.089 

12 1042 - 1100 1.459 3.957 

13 1056 - 1114 0.785 1.607 

14 1043 - 1101 0.589 0.861 

15 1642 - 1700 0.653 2.027 

16 1641 - 1659 0.585 2.296 

17 2243 - 2251 1.147 3.846 

18 0446 - 0504 1.079 4.346 

19 0442 - 0500 1. 208 2.509 

20 1042 - 1100 0.830 2.991 

Zi 1115 - 1133 0. 7Z1 2.538 

Z2 1111 - 1129 1. 27a 2.398 

23 1110 - 1128 1.135 4.344 

24 0510 - 0528 1.489 7.089 

25 1105 - 1123 1.208 6.737 

26 1103 - 1121 0.822 3.265 

2/ 1FO5 = 4123 1.456 2.897 

28 1102 - 1120 la55 4.869 

29 1043 - 1101 1.809 3.982 

30 1046 - 1104 1.369 5.709 

31 1043 - 1101 0.958 3.364 


82 


SUMMARY OF H3B AND H3S FOR JUNE, 1949 


Time 


1345 
1404 
0808 
1411 
1547 
1400 
0800 
0600 
0600 
1442 
1448 
1453 
1734 
P23 
1148 
1155 
2208 
1005 
1010 
1010 
1 
1134 
1145 
223 
1106 
1136 


(PCT) 


- 1403 
- 1422 
- 0826 
- Wee 
- 1605 
- 1418 
- 0818 
- 0618 
- 0618 
- 1500 
- 1506 
- 1511 
- 1752 
- 1a 
- 1206 
- 12a 
- 2226 
- 1023 
- 1028 
- 1028 
- Whee 
Bcy 
- i298 
- 1241 
- 14 
- 1154 
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HBB (ft.) 


ooo cocrr OrococooOr KR KB ke ke KF OND kK OOK CO 


741 
.@25 
.467 
. 753 
. 004 
.062 
-907 
. 205 
.005 
.045 
. 506 
« Sie. 
. 262 
865 
39 
878 
102 
.486 
12 / 
. 646 
465 
023 
745 
922 
719 
.861 


Has .(ft.) 


Pre WE OOWOWUF BS WWW WOU UNONF WHY ANE ND LEH 


.873 
- 187 
» 241 
-448 
641 
647 
-488 
. 100 
. 109 
862 
026 
O32 
979 
.973 
nD 
. 282 
. $11 
.374 
. 163 
.364 
2357 
672 
.146 
2519 
85/7 
«hoo 


SUMMARY OF H3B AND H3S FOR JULY, 1949 


Date Time (PCT) H3B (ft.) H3S (ft.) 
1 1948 - 2006 0.660 1.509 
zZ 1400 - 1418 0.470 1.207 
5 1355 - 1413 0.866 2 143 
4 1407 - 1425 0.798 2.695 
5 1355 - 1413 0.644 2.017 
6 1824 - 1842 0.915 2-328 
7 1223 - 1241 0.319 3.498 
9 0246 - 0304 0.860 3.923 

18 1352 - 1410 0.494 3.419 
19 1715 - 1733. 0.130 6.459 
20 1122 - 1140 1.062 272 
21 1125 - 1143 0.733 3.343 
22 1540 - 1558 1.098 5.164 
23 1545 - 1603 0.752 2.586 
24 1547 - 1605 0.635 2.392 
25 0949 - 1007 0.877 6.080 
26 1322 - 1340 0.730 4.686 
27 1457 - 1515 0.983 4.561 
28 1505 - 1523 1.049 4.595 
29 1506 - 1524 0.957 2.536 
30 1827 - 1845 0.932 2.655 
31 1230 - 1248 0.681 2.153 
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SUMMARY OF H3B AND H3S FOR AUGUST, 1949 


Date Time (PCT) H3B (ft.) H3S (ft.) 
1 1242 - 1300 0.876 3.214 
2 1457 - 1515 0.652 1.713 
3 1436 - 1454 0.733 1.474 
4 1530 - 1548 0.698 3.553 
5 1528 - 1546 1 e7 a. 556 
6 1504 - 1522 0.696 1.299 
Ul 1505 - 1523 1.089 6.107 
8 1503 - 1521 1.000 5.539 
9 1522 - 1540 0.960 2.860 

10 1521 - 1539 0.484 4.662 
m1 1520 - 1538 0.419 22S 
12 1529 - 1547 0.565 1.867 
13 1530 - 1548 0.769 3.641 
14 1528 - 1546 0.887 3.190 
5 1529 - 1547 0.524 0.876 
16 1541 - 1559 0.324 1.207 
7 1541 - 1559 0.574 1.555 
18 1541 - 1559 Oma L> 2.128 
19 0942 - 1000 0.627 1.162 
20 1211 - 1229 0.567 1.480 
21 1210 - 1228 0.240 1.972 
22 1208 - 1226 0.281 2.128 
23 1131 - 1149 0.288 1.586 
24 1130 - 1148 0 25T5 1.088 
25 1131 - 1149 0.500 2.943 
26 1129 - 1147 0.721 2.445 
27 1137 - 1155 0.824 Zr ole 
28 1139 - 1157 0.941 5.365 
29 1140 - 1158 0.590 2.158 
30 1201 - 1219 0.663 1.711 
31 1141 - 1159 0.506 2.891 
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